4 



DECUS NO. 


8-642 


TITLE 


AUTOCO - AUTOCORRELATION FOR POOR PEOPLE 
(WITHOUT EAE) 


AUTHOR 


Theodore J. Glattke 


COMPANY 


DATE 



SOURCE LANGUAGE 


Stanford University School of Medicine 
Stanford, California 

June 19, 1973 

XPAL 


ATTENTION 

This is a USER program. Other than requiring that it conform to submittal and review standards, 
no quality control has been imposed upon this program by DECUS. 

The DECUS Program Library is a clearing house only; it does not generate or test programs. No 
warranty, express or implied, is made by the contributor. Digital Equipment Computer Users 
Society or Digital Equipment Corporation as to the accuracy or functioning of the program or 
related material, and no responsibility is assumed by these parties in connection therewith. 







DECUS Program Library Write-up 


DECUS NO. 8-642 


AUTOCO (• t h EAF simulator) 
Overal l Prog ram De sign 


The program seeks to obtain an autocorrelation function on a string of 
data of 512 points or less by computing a Pearson product-moment 
correlation coefficient between elements of the firing and those same 
elements ’’delayed" with respect to themselves. This particular type of 
analysis may play an important role in the evaluation of random or 
quasi-random functions which may contain periodic components, for 
example, EEC activity coincident with repeated stimuli. (See the 
following: Lee, Y.W., Stat istical Theory of Communicatio n, New York: 

John Wiley and Sons, Inc., (i960) or Dern, H. , and Walsh, J.B., 

Analysis of complex waveforms, in Physical Techniques in Biological 
Research , Wm. L. Nastuk, editor, New York: Academic Press (1963) 

Volume 6.) 


The figure below provides an example of the effectiveness of this 
technique. The top tracing shows an example of a sample of quasi¬ 
random noise extending from 400 through 10,000 Hz, and an associated 
autocorrelation function based upon 200 points of that sample. 

The tracings on the bottom are another example of the noise mixed 
with a 500 Hz sinusoid at 0 dB S/N ratio. The autocorrelation function 
on the right clearly shows periodicity related to the fundamental 
component of the 500 Hz element, though its presence is not apparent 
in the input data. 



Figure 1. Examples of autocorrelation functions based upon 200 
data points with 100 microsecond bin-width. The magnitude of the 
autocorrelation function excursions is determined in part by S/N 
ratio, sample size (re: period of the component which is being sought), 
etc. 
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The data for which this program was designed are typically gathered 
on our 138E A/D converter. Data are input to the computer in a 10-bit 
signed format. Since our display is limited to a 10—bit unsigned 
format (it is a VC8/I), the input data are converted to unsigned 
numbers having a mean at lOOOg and a range from OOOOg through 177/g. 

It is in this form that we compute means, etc., of data - vi z •> 
unsigned positive numbers ranging from 0-1777 with a nominal mean 
at lOOOg. This is schematized in the figure below. 
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Figure 2. Representation of data suitable for input to the 
autocorrelator. 


It is not necessary for a user to have this particular type of system, 
only that he have his data in the format shown on the right side of 
the figure prior to use of the autocorrelator. 

The steps for the correlation are as follows. 


The computation' formula for the Pearson r is given as: 
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But since the X and Y distributions are based on the same data: 
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3. In the zero time delay condition, the expression reduces to: 
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Since working with 3 or 4 digits to the right ol the binary point 
is troublesome, the working formula used here further modifies the 
coefficient by adding a baseline of 1000c (to accomodate the display) 
after multiplying the coefficient by 1000 8 to arrive at ”r”values 
running from 1777g (for a +1.000) through lOOOg (for a 0.000) to 0001 
(for -1.000). There is an 0.2$ rounding error inherent in this 
but it should be satisfactory for most applications. 


Program flow 

This might be followed best by referral to the program listing. 
Generally, the procedure is as follows: 

1. Clear core and setup counters, etc. 

2. Read 4 blocks of data into locations 4000-4777 from unit 1. 

3. Find the smallest datum point and then subtract this (X sma ^^-1) 
from all data points to arrive at the smallest mean 

without affecting the range and without having any zeros. 

(This step saves much computation time.) 

4. Get SUMX „ 

5. Get (SUMX; 

6. Get N*SUM X 2 

7. Form denominator (which stays constant for this set of data) 

8. Get N*SUM XY (where Y is actually X "delayed” by 0,...,n bins) 

9. Form numerator (which changes with each ”delay' M 

10. Obtain correlation and store at location 5000,...,5777. 

11. Increment "delay" value and go back to 8 if not finished. 

12. If done, display data and correlation result. 

13. Search switch register. If bit 11=1, store r on unit 2, 
using next available 4 blocks. If bit 0=1, get next data from 
unit T and recycle through program. 


CRITICAL LOCATIONS 


V 


All user communication with program is via switch register, 
the following locations may be of interest to him. 


As such, 


1. Number of data blocks (at 128io words each): The program assumes 
that the user has 512 10 data points, or four (4) blocks of data 
on the source tape. If he has less than that, he may change the 
following locations to accomodate for the size of his field: 


Location Normal 

0011 0004 

0476 1204 

0453 2404 


Mnemonic 

# of blocks to be read/written 
Read function word, change last 
six bits to appropriate ft blks. 
Write function word, change Iasi 
six bits as above 


4 


3 
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2. Start of analysis — is assumed to be coincident with the first 
data point. You have the option of picking a new start point, 
e.g. , £-way down the sample. 

Location Normal Mnemonic 

0021 4000 Rpoint (1st point restoring value) 


If, for example, stimulus artifacts occupy the first few bins of 
your data, you may wish to set the starting point at 4100. 

The display resulting from this change in starting point will 
accomodate itself to it, and will locate the function over the 
data from which the function was derived. 

If you wish to change the starting point, keep the following 
in mind: 

3. The number of points may always be less than 512, but the 
Total number -- plus the starting point may never exceed 47771 

(or the end point of your data if you read in less than 4 blks) 
The 4777 is the address of the last location in which data are 
to be stored. If you exceed this, you'll include zeros in your 

analysis - a time consuming process. The counter is located as 

follows: 

Location Normal Mnemonic 

0023 7000 Rcount (2's complement of 

number of points in 
analysis) 

4. The starting block numbers for reading (location 0010) and 
writing (location 0012) may be changed freely, but the program 
assumes source data are on unit 1 and the autocorrelation 
functions will be stored on unit 2. Change the function word 
at 0453 (2404) to 1404 if you wish to store the functions on 
unit 1 (or any other by setting #404). 

. #* 





4 








AUTOCO page 5 
USE and ABUSE 


1. Loading and saving: You may load this program with a 1-pass 
binary loader. Core requirements are 0-1777, and starting address 
is 0200. 

2. After program is loaded, you may proceed as follows: 

A. Put source data on tape unit number 1. 

B. Put blank (marked) destination tape on unit number 2. 

C. Set new RPOINT and RCOUNT values, if desired. 

D. Load the starting address 0200 

E. and press START 

The program will run automatically until the function is completed 
for one-half of the total number of points in the sample, which 
corresponds to positive delay times re: the midpoint of the sample 
used for the analysis. In any non-random data, the continuation of 
the analysis would result in a mirror-image function corresponding 
to the"negative" delay times, and is unnecessary for most applications. 

The entire set of data read in will then be displayed, and the 
autocorrelation function will be superimposed upon the data, located 
on the x-axis at a point corresponding to the starting point of the 
analysis. 

F. If you wish to automatically store the data and proceed 
to obtain the next data , set the switch register to 4001. 

G. If you wish to store the data, but return to the display, 
set the switch register to 0001. You may start the next 
analysis by restarting the program at 0200. 

Using F above, you may put several samples of data on a source tape, 

start the first analysis, and go home - hopefully to find that 

several hundred functions have been completed by the morning. 

3. The range permitted is 0000-1777g for input data, but remember 
that the smaller the range, the shorter the computation time. 

In addition, there can be jno ZEROS as data points, for computing 
0 2 is very time-consuming! 

4. Our "means" are normally at 10008 because of the D/A system, but 
any non-zero mean less than 17778 should be satisfactory. The mean 
will not affect the autocorrelation output. 

5. Examples of the use of the program are attached. On the right in 
the figure, you may see five response means for stimuli presented to 
the auditory system at rates of 2000, 1000, 500, 250, and 125 pulses 
per second, respectively. 

The question we have regarding the responses is: "Is the cyclical 
activity shown between the vertical lines intersecting all the 
tracings related to the stimulus fundamental period?" 
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Figure 3. Examples of autocorrelation 
functions for brainstem following 
activity. 
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The autocorrelation functions corresponding to the data in the 
right column are shown in the left column. Normally, a relative 
or absolute maximum or minimum in the function relates to the 
period of some periodic event in the data; and zero axis crossings 
in the functions may also be interpretable in the same fashion. 

The arrows in each section of the autocorrelation functions 
show the anticipated locus of the period of the cyclical activity 
in the data corresponding to each function. As is apparent, 
one-to-one following activity was obvious for 250 and 125 pulses 
per second, but not for higher rates of stimulation. 

To familiarize himself with the use of the program, the user may 
prefer to analyze sinusoids or other repetitive stimuli in his 
laboratory with known temporal characteristics. 


Misc 


Autocorrelation functions are often displayed as indicated on the 
left below, viz., with "+" and delay times relative to a sample 

starting point. In a periodic signal, these should be mirror images, 
thus, in the present program, only the "+" delay results are computed 
and displayed. This is indicated on the right. 



Conventional plotting 



About time - the temporal sensitivity of the autocorrelation function 

output is identical to the sampling rate for the data which are submitted 
to analysis. If you are using a 100 microsecond/bin sample rate, then the 
points of the analysis will be equivalent to time delays of 100 microseconds. 
This is so because the delays are taken in one-bin intervals. If you 
plot your data on a conventional XY plotter and then plot the analysis 
results without changing the X-axis sensitivity, the time scale for the 
two displays will be identical. 

Computation time - this ranges from less than a minute to 17 minutes. 

It depends upon the number of points submitted to analysis and the 
absolute value of those points (squaring large numbers takes longer). 

You may wish to use only 32 or 64 points initially, until you have some 
facility with the program. 
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AUTOCORRELATION ROUTINE 
RECORD 1 


MAY >19 73 


T♦ GLATTKE 


OOfll 

3002 

0003 

0004 

0010 

0006 

DXL= 613 53 
DYL=6363 
*3010 

BEGBLK , 0003 

READ BLK NEXT (UNIT #1) 

0805 

001 1 

6004 

BLKS> 0004 

# BLKS TRANSFERRED 

0006 

0012 

0000 

WBLK, 0033 

WRITE BLK NEXT (UNIT #2) 

3007 

0013 

0000 

COKCNT, 3300 

# OF r VALUES (§ # points) 

0008 ■ 

0014 

.3000 

SIGN, 0003 

FOR SIGNED SUBTRACTION 

0009 

0015 

0000 

NEGATE, 0003 

SIGN FOR r 

0310 

03 1 6 

1000 

K 1 Id 0 0 , 1003 


0011 

031 7 

0777 

K 7 77, 0777 


0012 

0013 

3014 

\ 

0320 

4000 

*0023 

POINT, 4000 


00 1 5 

0021 

4000 

RPOINT, 4000 

1st DATA POINT (May be changed) 

0016 

0022 

7000 

COUNT, 7000 


0017 

0023 

7000 

RCOUNT, 7093 

2's comp of N (May be changed) 

0018 

0324 

5030 

CORR, 5090 


0019 

0025 

5000 

RCORR, 5060 

1st STORAGE LOCATION FOR r 

0020 

0026 

0000 

DPOINT, O000 

Incremented to cause "delay" 

0021 

0027 

0000 

CHECK, 0030 

Limit of delay 

0022 

0330 

0000 

DELAY, 3090 

Increment and added to DPOINT 

0fP23 

0024 

0031 

7000 

i 

RESET, NOP 


0025 

0332 

7300 

CLA CLL 

Reset pointer and counter 

0026 

0033 

1321 

TAD RPOINT 

0027 

0334 

3620 

DC A PO INI- 


0028 

003 5 

1 623 

TAD RCOUNT 


6029 

0336 

3022 

DCA COUNT 


0030 

003 7 

3151 

DCA XAXIS 


6031 

0040 

5431 

JMP I RESET 


0032 

3333 

0641 

7000 

SETCOR, NOP 


0034 

00 42 

7360 

CLA CLL 

Reset Storage for r 

0035 

0043 

1025 

TAD RCORR 


0036 

00 4 4 

3024 

DCA CORR 


0037 

0645 

5441 

JMP I SETCOR 


3033 

0339 

0046 

7000 

CORE, NOP 


0340 

004 7 

41 53 

JMS XSET 

Clears 20 OO 3 locations beginnin 

03 41 

03 50 

3420 

DCA I POINT 

at 4000, and going through 5777 

0042 

00 51 

3424 

DCA I CORR 

00 43 

00 52 

2023 

ISZ POINT 


00 44 

03 53 

26 24 

ISZ CORR 


0345 

00 54 

2022 

ISZ COUNT 


0046 

60 55 

56 53 

JMP .- 5 


00 4 7 

03 56 

4331 

JMS RESET 


0043 

0O5 7 

3030 

DCA DELAY 


00 49 

00 6O 

5446 

JMP I CORE 


00 50 

00 51 

00 61 

7000 

LIMIT, NOP 

Sets "delay” = RPOINT+N at its 

00 52 

•30 62 

1323 

TAD RCOUNT 

Caution Ill 

00 53 

00 63 

70 4 1 

CIA 

This can never be greater than 
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AUTOCOFRELA TI ON ROUT INE 
RECORD 1 


MAY#1973 


T. GLATTKE 


00 54 
0055 
0956 
0057 
00 53 
00 59 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
,0067 

^p)063 
00 69 
00 70 
00 71 
09 72 
0073 
0074 
0075 
00 76 
0077 
09 78 
00 79 
00 80 
0081 
0082 
0083 
0084 
0085 

fP3087 
9088 
0089 
0090 
0091 
0092 
0093 
0094 
009 5 
0096 
039 7 
0098 
0099 
9100 
3131 
0102 
3133 
0104 
01 3 5 
010 6 
0107 


0964 
0065 
0066 
336 7 


1021 TAD RPOINT 
7041 CIA 
3027. DCA CHECK 
5461 JMP I LIMIT 


00 70 7000 

0071 7300 

0072 1023 

00 73 7041 

0874 71 10 

00 75 70 41 

0076 3013 

0077 5470 


HALF# NOP 
CLA CLL 
TAD KCOUNT 
CIA 

CLL RAR 
CIA 

DCA CORCNT 
JMP I HALF 


Sets r counter = \ number of points. 
To get full function, insert NOPs 
in 73,74,75. 

In this mode, only the "+" delay 
times are used. 


8109 8000 

0101 0008 
0102 0000 
0103 0000' 

3104 0030 

0105 0000 

0106 0000 
0107 0003 

0110 0000 
0111 0003 
0112 0000 
0113 0003 

0114 0000 

0115 0080 

0116 0000 
0117 0000 

0120 0000 
0121 0080 
0122 0300 
0123 0800 

0124 0000 

0125 0000 

0126 0000 
0127 0000 

0130 0903 

0131 0000 

0132 0303 

0133 0000 


SUM 1 #0330^ 
SUM?#0000 
SUM3# 0000 1 
SUM4# 0008 J 
OR 1#0000S 
0R2# 0000 I 


0R3 # 0000 
0R4# 3000 
XI #0003 
X 2 # 0 0 (5 8 
X3# 0980 


X4#0000 
Y 1 #0088 
Y2# 0080 


Y 3#8000 I 
Y 4 # 0003 0/ 
SMX 21 # 0000 N 
SMX 22#0000 L 
SMX 23 # 0000 T 
SMX 24 # 0080 J 
DENOM!# 0000 S 
DEN 0 M 2 # 0300 t 
DEN 0 M 3#0003 C 
DENOM 4#0000 / 


XY1 #0030 
XY2# 0300 
XY3# 0000 I 
XY 4# 0000 J 


Registers for 



Accumulators for multiply and 
subtraction routines 


Registers for 


Registers for 


Registers for 


K (it*) 

Kt 


013 4 

0472 

READ# 

GET 

0135 

0534 

SUMX# 

GETX 

3136 

0 690 

MUL# MPY 

0137 

1030 

SUMX2# 

GETX2 

3 1 40 

1 1 22 

DENOM# 

DIVIS 

0141 

0 716 

SUBTR# 

MINUS 

0142 

1 480 

XPROD# 

XYPEOD 

91 43 

1 680 

NUMF.R # 

UI VI D 

0 1 44 

0088 

TALLY# 

9990 

01 4 5 

0090 

TF LAG# 

0000 


Indirect address tags and 
misc stuff 
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AUTOCORRELATION ROUTINE MAY*1973 T. GLATTKE 

RECORD 1 


0103 

01 46 

0002 

I NCR* 0002 

0109 

01 47 

1200 

RW* RWTAPE 

0110 

0 1 50 

4000. 

KPOINT* 4000 

0111 

0151 

0000 

XAXIS* 0000 

0112 

0152 

•3000 

KEEP* 0030 

0113 

0153 

7000 

XSET* NOP 

0114 

0154 

1 161 

TAD K7 

01 1 5 

01 55 

3022 

DCA COUNT 

01 1 6 

0156 

1162 

TAD K2 

01 1 7 

0157 

3020 

DCA POINT 

01 18 

01 60 

5553 

JMP I XSET 

01 19 

0161 

7000 

K 7 * 70013 

0120 

0121 

0162 

4000 

K2* 4000 
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MUT0C0RK ELA T10INI ROUT I N F. 
RECORD 2 


MAY*1973 


To GLATTKE 


000 1 
0002 

0209 

7300 

*0200 

START, CLA CLL 

0303 

0201 

40 41 

JMS 

SETCOR 

0304 

0202 

4046 

JMS 

CORE 

0305 

0203 

40 70 

JMS 

HALF 

903 6 

020 4 

4534 

JMS 

I READ 

033 7 

023 5 

4031 

JMS 

RESET 

0038 

0236 

4041 

JMS 

SETCOR 

3009 

020 7 

4153 

JMS 

XSET 

0010 

0210 

1 420 

TAD 

I POINT 

001 1 

021 1 

70 41 

CIA 


09 1 2 

0212 

3927 

DCA 

CHECK 

0013 

0213 

1 420 

TAD 

I POINT 

001 4 

0214 

1027 

TAD 

CHECK 

0315 

021 5 

7500 

SMA 


001 6 

0216 

5223 

JMP 

. + 5 

9017 

021 7 

7300 

CLA 

CLL 

0018 

0220 

1 429 

TAD 

I POINT 

0019 

0221 

7041 

CIA 


0020 

0222 

3027 

DCA 

CHECK 

0021 

0223 

7300 

CLA 

CLL 

0022 

0224 

2020 

ISZ 

POINT 

0023 

9225 

2022 

ISZ 

COUNT 

0324 

3226 

5213 

JMP 

.-13 

0325 

0227 

41 53 

JMS 

XSET 

0026 

0230 

1027 

TAD 

CHECK 

0027 

0231 

7001 

I AC 


3928 

9232 

3027 

DCA 

CHECK 

0329 

0233 

1 420 

TAD 

I POINT 

3033 

3234 

1027 

TAD 

CHECK 

3031 

0235 

3429 

DCA 

I POINT 

0032 

0236 

2020 

ISZ 

POINT 

0033 

0237 

2022 

ISZ 

COUNT 

0034 

0240 

5233 

JMP 

.-5 

0935 

0241 

4031 

JMS 

RESET 

0936 

0242 

4061 

JMS 

LIMIT 

93 37 

0038 

0243 

4535 

JMS 

I SUMX 

3339 

0040 

0244 

4536 

JMS 

I MUL 

0041 

0245 

1 104 

TAD 

OR 1 

0342 

0246 

3100 

DCA 

SUM i 

9043 

0247 

1 1 35 

TAD 

OR 2 

0944 

•32 59 

3101 

DCA 

SUM2 

09 45 

0251 

1 106 

TAD 

OR 3 

0046 

0252 

3102 

DCA 

SUM 3 

09 4 7 

3253 

1197 

TAD 

0R4 

3048 

92 54 

3103 

DCA 

SUM 4 

0949 

39 59 

0255 

7300 

CLA 

CLL 

00 51 

m 52 

0256 

453 7 

JMS 

I SUMX2 


A,) 53 


Restore r pointer 

Clear core from 4000-5777 

Set r counter to \ N (optional) 

Get DATA from tape unit #1 

Set pointers and counters 

Restore r pointer 

Set counter for 512 ^q points 


Find smallest datum point (X s ) 


NOTE NOTE NOTE-THERE CAN BE 

NO 0000s IN YOUR SAMPLE!!!! 


Subtract (X s -1) from all data 
so that the mean is reduced to 
smallest possible value, but 
the variance is not disturbed. 


Puts into SUM 1 ... SUM3 

Puts/^ hV' into SUM1 ... SUM3 


Puts into SUMX21 ... SUMX23 
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AUTOCORRELATION ROUTINE 
RECORD 2 


MAY j> 1973 


T. GLATTKE 


00 54 
0055 

0257 

4540 

JMS 

I DENOM 

0056 

0057 

0263 

73.00 

CLA 

CLL 

0058 

0261 

4031 

JMS 

RESET 

0059 

0060 

0262 

3030 

DCA 

DELAY 

0061 

0062 

0263 

4542 

CROSPD, JMS I 

0063 
00 64 

0264 

7300 

CLA 

CLL , 

0065 

0366 

0265 

4543 

JMS 

I NUMER 

0067 
00 68 

0266 

3144 

DCA 

TALLY 

03 69 
00 70 

0267 

7303 

CLA 

CLL 

00 71 

02 73 

1 104 

DOCOj TAD ORl 

00 72 

0271 

3110 

DCA 

XI 

33 73 

0272 

110 5 

TAD 

QR2 

30 74 

0273 

31 1 1 

DCA 

X2 

03 75 

0274 

1 106 

TAD 

0R3 

0076 

0275 

3112 

DCA 

X3 

0377 

0276 

1 137 

TAD 

0R4 

30 78 

0277 

3113 

DCA 

X 4 

0-379 

0303 

1 124 

TAD 

DENOM1 

0383 

0301 

31 1 4 

DCA 

Y 1 

008 1 

033? 

1 125 

TAD 

DEN0M2 

0082 

0303 

31 1 5 

DCA 

Y 2 

0033 

3 3 0 4 

1 1 26 

TAD 

DEN0M3 

0384 

033 5 

31 1 6 

DCA 

Y3 

0385 

033 6 

1 127 

TAD 

DENOM4 

3086 

033 7 

3117 

DCA 

Y 4 

0387 

0310 

4541 

JMS 

I SUBTft 

0088 

031 1 

1014 

TAD 

SIGN < 

0389 

3312 

7650 

SNA 

CLA 

3090 

0313 

5321 

JMP 

• + 6 

0091 

0314 

7300 

CLA 

CLL 

0392 

3315 

1 1 44 

TAD 

TALLY 

0093 

0316 

7001 

I AC 


039 4 

331 7 

3144 

DCA 

TALLY 

0395 

0320 

5270 

JMP 

DOCO 

00 9 6 

0321 

7300 

CLA 

CLL 

039 7 

0322 

1015 

TAD 

NEGATE 

0398 

0323 

7440 

SZA 


0099 

0324 

5332 

JMP 

« + 6 

0100 

0325 

7300 

CLA 

CLL 

0101 

0326 

1 1 44 

TAD 

TALLY 

010? 

0327' 

1017 

TAD 

K777 

0103 

3 333 

3 1 4 4 

DCA 

TALLY 

0104 

0331 

5337 

JMP 

• + 6 

3105 

0332 

7330 

CLA 

CLL 

010 6 

0 333 

1 1 44 

TAD 

TALLY 

3 1 0 7 

0334 

7041 

CIA 



Puts 


N £(.**)•(**) 


into DENOM1 
thru DEN0M3 




into XY1 


XY3 


Forms 


f|4£*y * in OR1 


thru OR4 


Division is by repeated subtraction 


Done? 

Yes 


Get sign 

Form + correlation value 


form (-) correlation value 
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AUTOCORRELATION ROUTINE 

RECORD 2 


MAY,1973 


T. GL.ATTKE 


0108 

0335 

1016 

TAD 

K 1 000 

0109 

0336 

3144 

DCA 

TALLY 

01 10 

033 7 

1 1 44 

TAD 

TALLY 

01 1 1 

0340 

3424 

DCA 

I CORR 

01 12 

0341 

2024 

ISZ 

CORR 

01 13 

0342 

2030 

ISZ 

DELAY 

01 14 

03 43 

3015 

DCA 

NEGATE 

011 5 

0344 

2013 

ISZ 

CORCNT 

01 1 6 

0345 

5263 

JMP 

CROSPD 

01 17 

0346 

7303 

CLA 

CLL 

01 18 

0347 

3145 

DCA 

TFLAG 

0119 

0350 

5751 

JMP 

I NEXT 

0120 

0121 

0351 

0400 

NEXT 

0400 
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AUTOCORRELATION' ROUTINE MAY, 1973 
RECORD 3 


000 1 



*0430 

0002 

0400 

7300 

DISPLYj cla cll 

0003 

0481 

430 7 

JMS 

RESETD 

609 4 

0402 

1 151 

TAD 

XAXIS A 

0005 

0403 

6953 

DXL 


(3006 

0404 

7200 

CLA 


900 7 

0405 

1 423 

TAD 

I POINT 1 

0008 

0406 

63 63 

DYL 


0909 

8 40 7 

720-3 

CLA 


8010 

0410 

2020 

ISZ 

POINT 1 

08! 1 

041 1 

1 1 51 

TAD 

XAXIS 

00 1 2 

0412 

1 1 46 

TAD 

i ncr r 

0013 

0413 

31 51 

DCA 

XAXIS 1 

801 4 

8)41 4 

2022 

ISZ 

COUNT 

001 5 

0415 

5232 

JMP 

♦ - 1 3 i 

0816 

@416 

7300 

CLA 

CLL *** 

0817 

041 7 

4317 

JMS 

RESETC 

0018 

0420 

1 151 

TAD 

XAXIS \ 

0919 

0421 

60 53 

DXL 

\ 

0020 

04 ?? 

7280 

CLA 

\ 

0821 

0423 

1 420 

TAD 

I POINT \ 

0022 

0424 

70 1 0 

RAR 

v 

0923 

0425 

6063 

DYL 

r 

0024 

0426 

7206 

CLA 

1 

•3025 

0427 

2820 

ISZ 

POINT 1 

0026 

0433 

1 1 51 

TAD 

XAXIS 

0027 

0431 

1 1 46 

TAD 

I NCR I 

0823 

0432 

31 51 

DCA 

XAXIS J 

0029 

0433 

2022 

ISZ 

COUNT +S 

0930 

0434 

5220 

JMP 

«-l 4 

9931 

0435 

7333 

CLA 

CLL 

0032 

0436 

7604 

LAS 


0033 

043 7 

70 1 0 

RAR 


0334 

8 448 

7620 

SNL 

CLA 

8835 

8441 

5200 

JMP 

DISPLY 

0936 





0837 

0442 

7300 

CLA 

CLL 

0838 

6443 

1 1 45 

TAD 

TFLAG 

0039 

0444 

7440 

SZA 


00 40 

0 445 

5230 

JMP 

DISPLY 

3341 

0446 

7001 

I AC 


30 42 

0447 

3145 

DCA 

TFLAG 

00 43 

0450 

1812 

TAD 

WBLK 

00 44 

0451 

3254 

DCA 

OUT 

•3045 

8 4 52 

4547 

JMS 

I RW 

0046 

0453 

2404 

2 40 4 

0047 

0454 

0000 

OUT > 

0 

9948 

0455 

5900 

5900 

i 

3949 

0 4 56 

8457 

. +1 


89 50 

0457 

1812 

TAD 

WBLK 

09 5 1 

3 460 

10 11 

TAD 

BLKS 

33 52 

0461 

3812 

DCA 

WBLK 

9353 

0 462 

7300 

CLA 

CLL 


T. GLATTKE 


Displays data and r function 
on VC8/I. 

Note, XAXIS is incremented by 
2 to accomodate full-width display. 

User may wish to modify to put 
data on left and e.jg., r on right. 

Displays contents of 4000-4777. 
(Data) 


Displays contents of 5000 and up, 
to a limit of 5000 + N;, and 
locates the display over the data 
at point where analysis is initiated. 


Bit 11 = 1? 



No, continue with display 
Yes, have data been stored? 


No, raise tape flag 

and put contents of 5000-5777 in 

next available block on Unit #2. 
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AUTOCORRELATION ROUTINE 
RECORD 3 


MAY* 1,9 73 


T, 6LATTKF 


00 54 

3 463 

763 4 

LAS 

00 55 

3464 

730 4 

RAL 

00 5 A 

346 5 

7620 

SNL CLA 

00 5 7 

0 466 

5230 

JMP D1SPLY 

00 58 

046 7 

56 70 

JMP I XSTART 

0059 

0 4 70 

0200 

XSTART * 3200 

00 60 
00 61 

0471 

7000 

NOP 

00 62 

0472 

7000 

GET* NOP 

0363 

04 73 

1013 

TAD BEGBLK \ 

0064 

0474 

3277 

DCA INB 

0065 

0475 

4547 

JMS I RW 1 

0066 

04 76' 

1 204 

1204 \ 

00 6 7 

04 77 

0000 

INB*3 f 

0068 

0503 

4003 

4009 

0069 

0 501 

0502 

♦ + 1 1 

03 70 

0 532 

1010 

TAD BEGBLK 

0071 

0503 

1011 

TAD BLKS 1 

m 72 

0504 

3010 

DCA BEGBLK 1 

00 73 

050 5 

7303 

CLA CLL 1 

00 74 
0075 

0536 

5672 

JMP I GET 

0076 

0 53 7 

7003 

RESETD* NOP 

mu 

0510 

7300 

CLA CLL 

00 78 

051 1 

1 1 53 

TAD KPOINT 

00 79 

0512 

3020 

DCA POINT 

0083 

0513 

1333 

TAD DXNT 

0081 

351 4 

3022 

DCA COUNT 

0082 

051 5 

3151 

DCA XAXIS 

0083 

0516 

573 7 

JMP I RESETD 

0084 

051 7 

7303 

RESETC* NOP 

008 5 

0520 

7330 

CLA CLL 

003 6 

0521 

1025 

TAD RCORR 

0087 

0522 

3023 

DCA POINT 

0088 

0523 

43 70 

JMS HALF 

0089 

0524 

1 0 1 3 

TAD CORCNT 

0090 

3 52 5 

3022 

DCA COUNT 

0091 

0 526 

1021 

TAD RPOINT 

0392 

0527 

70 3 4 

RAL 

0093 

3 530 

3151 

DCA XAXIS 

039 4 

3531 

5717 

JMP I RESETC 

0095 

0532 

4000 

K4* 4030 

309 6 
009 7 

0533 

7303 

DXNT* 7030 


After data transfer, check bit 0 
Bit 0=1? 

No, go to display 
Yes, go to program start 

Entry to tape read routine, 
gets next 4 blocks from Dectape #1 
and puts data into loc 4000-4777 
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AUTOCORRELATION ROUTINE 
RECORD 4 


MAY#1973 


T. GLATTKE 


1 



mdi 

0534 

7380 

GETX 

, NOP 

0002 

0535 

7330 

CL A 

CLL 

0003 

0536 

3100 

DCA 

SUM 1 

0004 

0537 

3101 

DCA 

SUM2 

0005 

0540 

3132 

DCA 

SUM3 

000 6 

0541 

3103 

DCA 

SUM4 

0007 

0542 

4331 

JMS 

RESET 

0008 

0 543 

1 420 

TAD 

I POINT 

0009 

0544 

1 100 

TAD 

SUM 1 

0010 

0545 

3100 

DCA 

SUM 1 

001 1 

0546 

730 4 

RAL 


001 2 

0 54 7 

1101 

TAD 

SUM2 

0013 

0 550 

3101 

DCA 

SUM2 

0014 

0551 

2020 

ISZ 

POINT 

0015 

0552 

2022 

ISZ 

COUNT 

0016 

0553 

5343 

JMP 

.-10 

001 7 

0554 

4031 

JMS 

RESET 

0018 

0555 

1 1 00 

TAD 

SUM 1 

0019 

0556 

3119 

DCA 

XI 

0020 

0557 

1101 

TAD 

SUM2 

0021 

8560 

31 1 1 

DCA 

X2 

0022 

0561 

1 100 

TAD 

SUM 1 

0023 

0562 

31 1 4 

DCA 

Y1 

0024 

8563 

1101 

TAD 

SUMP 

0025 

0564 

3115 

DCA 

Y2 

0026 

3565 

3112 

DCA 

X3 

0027 

0566 

3113 

DCA 

X4 

0028 

0567 

3116 

DCA 

Y 3 

0029 

0570 

31 1 7 

DCA 

Y 4 

0030 

0571 

5734 

JMP 

I GETX 

0031 



*0600 

0332 





0033 

0600 

7000 

MPYj 

. NOP 

0034 

0601 

7303 

CLA 

CLL 

003'5 

0602 

3104 

DCA 

OR 1 

0036 

0603 

3105 

DCA 

0R2 

0037 

0604 

3106 

DCA 

OR3 

0038 

0605 

3107 

DCA 

OR 4 

0039 

0 60 6 

1114 

TAD 

Y1 

0340 

063 7 

73 41 

CIA 


0041 

0610 

3114 

DCA 

Y1 

0042 

061 1 

7004 

RAL 


0043 

0612 

3152 

DCA 

KEEP 

0044 

0613 

1115 

TAD 

Y2 

0045 

3 61 4 

73 40 

CMA 


0046 

061 5 

1 1 52 

TAD 

KEEP 

004 7 

0616 

3115 

DCA 

Y2 

00 48 

-361 7 

7004 

RAL 


0049 

0620 

31 52 

DCA 

KEEP 

00 50 

0 621 

1116 

TAD 

Y 3 

03 51 

0622 

7 0 43 

CMA 


30 52 

0 623 

1 1 52 

TAD 

KEEP 

0053 

0624 

3116 

DCA 

Y3 


Gets U and puts into 
X and Y registers for squaring 


General multiplication routine for 
48 bit unsigned numbers. 

Clear MUL registers 


Form 2's complement of multiplier 
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AUTOCORRELATI ON ROUT I NF. 
RECORD 4 


MAY,1973 


T. GLATTKE 


0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 
0063 
0064 
0065 
0966 
0067 
0068 
00 69 
00 70 
00 71 
00 72 
0973 
00 7 4 
00 75 
0076 
0077 
30 73 
00 79 
0030 
0081 
0082 
0083 
0084 
0085 
0036 
0387 
0088 
0089 
0090 
0091 
0092 
0093 
9394 
0095 
909 6 
0097 
0098 
0099 
0100 
01 3 1 
0102 
0103 
0134 . 
0105 
0106 
0107 . 


0625 

730 4 

0626 

31 52 

0 627 

1117 

0 630 

7040 

0631 

1 1 52 

0632 

3117 

0633 

7300 

0634 

1110 

0635 

1104 

0636 

3104 

3 63 7 

700 4 

0 643 

1 105 

0641 

3195 

06 42 

700 4 

0643 

1 106 

0644 

3106 

0645 

7004 

0646 

1107 

0647 

3107 

0 650 

7300 

3651 

1111 

0652 

1105 

06 53 

3105 

0654 

700 4 

0655 

1 106 

0656 

3106 

@657 

7034 

0 660 

1107 

0661 

3107 

0 662 

7300 

06 63 

1112 

0 664 

1136 

0665 

3106 

0 666 

7004 

0667 

1 107 

0 670 

3107 

0671 

7390 

0672 

1113 

0673 

1107 

0 6 7 4 

310 7 

0675 

7390 

0676 

7001 

0677 

1114 

0 700 

3114 

0 70 1 

700 4 

0702 

1115 

0 733 

3115 

@70 4 

7304 

0 73 5 

1116 

0 73 6 

3116 

3 70 7 

7004 

0 710 

1117 

0 71 1 

3117 

0712 

1117 


RAL 

DCA KEEP 
TAD Y4 
CMA 

TAD KEEP 
DCA Y 4 
DO* CLA CLL 
TAD XI 
TAD OR 1 
DCA OR 1 
RAL 

TAD 0R2 
DCA 0R2 
RAL 

TAD 0R3 
DCA 0R3 
RAL 

TAD 0R4 
DCA OR4 
CLA CLL 
TAD X2 
TAD 0R2 
DCA 0R2 
RAL 

TAD 0R3 
DCA 0R3 
RAL 

TAD 0R4 
DC A 0K4 
CLA CLL 
TAD X3 
TAD 0R3 
DCA 0R3 
RAL 

TAD OR4 
DCA 0R4 
CLA CLL 
TAD X4 
TAD OR4 
DCA OR4 
CLA CLL 
I AC 

TAD Y 1 
DCA Y1 
RAL 

TAD Y2 
DCA Y2 
RAL 

TAD Y3 
DCA Y3 
RAL 

TAD Y4 
DCA Y 4 
TAD Y4 


Increment OR registers with value 
of mutliplicand until 2’s complement 
of multiplier is = 0. 

The product of this is in 0R1 ... 0R3 


17 





AUTOC0 RR ELATION ROUTINE 
RECORD 4 


MAY > 1973 


T. GLATTKE 


t 


0108 

3713 

74 53 

SNA 


3109 

0714 

5603 

JMP 

I MPY 

31 10 
0111 

071 5 

5233 

JMP 

DO 

01 12 

0716 

7030 

MINUS, NOP 

01 13 

0 717 

3014 

DC A 

SI GN 

0114 

0 720 

1114 

TAD 

Yi 

0115 

0721 

70 41 

CIA 


01 1 6 

0722 

1113 

TAD 

XI 

0117 

0723 

3104 

DCA 

OR! 

0118 

3 724 

7904 

RAL 


01 19 

0725 

3152 

DCA 

KEEP 

0120 

3 726 

1115 

TAD 

Y2 

3121 

0727 

7049 

CMA 


0122 

3 730 

1 1 1 1 

TAD 

X2 

0123 

0 731 

1 1 52 

TAD 

KEEP 

0124 

3732 

3105 

DCA 

0R2 

0125 

0 733 

7004 

RAL 


0126 

0 734 

3152 

DCA 

KEEP 

0127 

9 735 

1116 

TAD 

Y3 

0128 

0 736 

70 40 

CMA 


0129 

0737 

1112 

TAD 

X3 

0139 

3 749 

1 1 52 

TAD 

KEEP 

3131 

0 741 

3106 

DCA 

0R3 

0132 

0 742 

7004 

RAL 


0133 

0 743 

31 52 

DCA 

KEEP 

0134 

9 744 

1117 

TAD 

Y 4 

0135 

0 745 

70 40 

CMA 


3136 

0 746 

1113 

TAD 

X 4 

0137 

0 74 7 

1 1 52 

TAD 

KEEP 

0138 

0753 

319 7 

DCA 

OR 4 

0139 

0751 

7094 

RAL 


0143 

0 752 

3014 

DCA 

SIGN 

01 41 

0142. 

0 753 

5716 

JMP 

I MINUS 


Form 2's complement of subtrahend 
and add it to minuend. Both numbers 
are unsigned and may be up to 48 bits 



SIGN= 1 if result is positive 
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AUTOCORRelation routine 
RECORD 5 


MAY*1973 


T. GLATTKE 


000 1 



♦ 1 030 


0002 

1000 

7300 

GETX2, NOP 

1 

^ *J O 

1001 

4331 

JMS RESET 


0004 

1 002 

31 13 

DCA XI 


000 5 

1033 

31 1 1 

DC A X2 


000 6 

133 4 

3112 

DCA X3 


0007 

1335 

3113 

DCA X4 

Clear multiplier, multiplicand, 

0008 

103 6 

31 1 4 

DCA Y1 

and accumulator registers for the 

0039 

100 7 

3115 

DCA Y2 

XY cross-products 

00 1 0 

10 10 

3116 

DCA Y3 


001 1 

1011 

31 1 7 

DCA Y4 


0312 

1012 

3120 

DCA SMX21 


0013 

1013 

3121 

DCA SMX22 


00 1 4 

1014 

3122 

DCA SMX23 


m i 5 

1015 

3123 

DCA SMX24 


0016 
03 1 7 

1016 

1 420 

PNTR, TAD I 

POINT 

0318 

1017 

31 10 

DCA XI 

> 

3319 

1020 

1 420 

TAD I POINT 

Increment pointer, squaring the 

0023 

1021 

31 1 4 

DCA Y1 

contents of the location indicated 

0321 

1022 

4536 

JMS I MUL 

by it until all relevant points 

0-322 

1023 

7303 

CLA CLL 

have been squared and summed in 

0323 

1024 

1 104 

TAD OR 1 

SMX21 ... SMX23 

0324 

1025 

1 1 20 

TAD SMX21 


0325 

1 026 

3120 

DCA SMX21 


0326 

1 027 

733 4 

RAL 


3027 

1033 

1121 

TAD SMX22 

• 

0323 

1331 

3121 

DCA SMX22 

• 

0229 

1032 

73£M 

RAL 


0030 

1033 

1 122 

TAD SMX23 


0331 

1034 

3122 

DCA SMX23 


0032 

1035 

7304 

RAL 


0333 

1036 

1 1 23 

TAD SMX24 


0034 

1037 

3123 

DCA SMX24 


0035 

1043 

7300 

CLA CLL 


0036 

1(341 

1105 

TAD 0R2 


003 7 

1342 

112! 

TAD SMX22 


0338 

1343 

3121 

DCA SMX22 


0039 

1044 

700 4 

RAL 


00 40 

1045 

! 1 22 

TAD SMX23 


00 41 

1346 

3122 

DCA SMX23 


00 42 

1047 

7004 

RAL 


00 43 

13 53 

1 123 

TAD SMX24 


0044 

1351 

3123 

DCA SMX24 


00 45 

10 52 

7330 

CLA CLL 


00 46 

1353 

1 106 

TAD 0R3 


034 7 

13 54 

1 1 22 

TAD SMX23 


0043 

1 0 55 

3122 

DCA SMX23 


00 49 

1 3 56 

70 0 4 

RAL 


•-0 c; o 

13 57 

1 1 23 

TAD SMX24 


•30 51 

1 3 63 

31 23 

RCA SMX24 


03 52 

1 361 

7303 

CLA CLL 


00 S3 

10 62 

1107 

TAD OR4 



19 




AUTOCORRELATION ROUTINE 
RECORD 5 


MAY > 19 73 


T. GLATTKE 


^ 00 54 
'0055 

10 63 

1 123 

TAD 

SMX24 

10 64 

3123 

DCA 

SMX24 

0056 

1065 

2020 

ISZ 

POINT 

0057 

1066 

20 22 

ISZ 

COUNT 

0058 

1067 

5216 

JMP 

PNTR 

0059 

1070 

4031 

JMS 

RESET 

00 60 

1071 

l 1 20 

TAD 

SMX21 

0061 

1072 

31 10 

DCA 

XI 

0062 

10 73 

1 121 

TAD 

SMX22 

0063 

1074 

31 1 1 

DCA 

X2 

0064 

1075 

1 122 

TAD 

SMX23 

0065 

1076 

31 12 

DCA 

X3 

0066 

10 77 

1 123 

TAD 

SMX2 4 

0967 

1 100 

31 13 

DCA 

X 4 

00 68 

1101 

1023 

TAD 

RCOUNT 

0069 

1 102 

70 41 

CIA 


00 70 

1 103 

3114 

DCA 

Y1 

00 71 

1 104 

31 1 5 

DCA 

Y2 

00 72 

1 105 

31 1 6 

DCA 

Y3 

00 73 

1106 

3117 

DCA 

Y 4 

00 74 

1107 

4536 

JMS 

I MUL 

00 75 

1110 

7300 

CLA 

CLL 

03 76 

1111 

1 104 

TAD 

OR 1 

00 77 

1112 

3120 

DCA 

SMX21 

00 78 

1113 

1 105 

TAD 

0R2 

0979 

1114 

3121 

DCA 

SMX22 

0080 

1115 

1 196 

TAD 

0R3 

9381 

1116 

3122 

DCA 

SMX23 

998? 

1117 

1197 

TAD 

OR 4 

9083 

1 120 

3123 

DCA 

SMX24 

0084 
308 5 

1 121 

5600 

JMP 

I GETX2 

0086 

1122 

7000 

DIVIS* NOP 

3087 

1 1 23 

7300 

CLA 

CLL 

0088 

1 124 

1 1 20 

TAD 

SMX21 

0089 

1 125 

31 10 

DCA 

XI 

0090 

1 1 26 

1121 

TAD 

SMX22 

0091 

1 127 

31 1 1 

DCA 

X2 

0092 

1 130 

1 1 22 

TAD 

SMX23 

0993 

1131 

3112 

DCA 

X3 

0094 

1 132 

1 123 

TAD 

SMX24 

9095 

1 133 

31 1 3 

DCA 

X 4 

0096 

1 134 

1 1 00 

TAD 

SUM 1 

0997 

1 135 

3114 

DCA 

Y 1 

0098 

3 1 36 

110 1 

TAD 

SUM 2 

0099 

113 7 

3115 

DCA 

Y? 

0100 

1 l 40 

1 1 02 

TAD 

SUM3 

919 1 

1 141 

3116 

DCA 

Y3 

0132 

1 1 42 

1 103 

TAD 

SUM 4 

0103 

1 1 43 

311"? 

DCA 

Y 4 

0194 

1 l 44 

4541 

JMS 

I SUBTR 

U I ~ 3 

01 0 6 

1 1 45 

1 1 0 4 

TAD 

OR 1 

010 7 

1 1 46 

3124 

DCA 

DENOM1 


Done with all points? 

No 

Yes, now put the sum of the 
squared scores into the multiplicand 


and put N into multiplier 



and get 




in the SMX2 registers 


Divisor is formed by subtracting 
from and the 


result is . stored in DENOM registers 
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AUTOCORRELATION ROUTINE 
RECORD 5 


MAY,1973 


T. 6LATTKE 


0108 

1 1 47 

1105 

TAD 

0R2 

3109 

1 1 53 

3125 

DC A 

DEN0M2 

3110 

1151 

1 106 

TAD 

0R3 

01 1 1 

1 1 52 

3126 

DCA 

DEN0.M3 

0112 

1 1 53 

1137 

TAD 

0R4 

0113 

1 1 54 

3127 

DCA 

DEN0M4 

01 1 4 
01 1 5 

1 1 55 

5722 

JMP 

I DIVIS 
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AUTOCORRELATION ROUTINE MAY,1973 
RECORD 6 


T. GLATTKE 


000 1 
0002 

1 200 

7432 

*1200 

RWTAPE, HLT 

0003 

1201 

7330 

CLA CLL 

0004 

1202 

1 600 

TAD I RWTAPE 

0005 

1203 

2200 

ISZ RWTAPE 

000 6 

1204 

3325 

DCA RWREGA 

0007 

1235 

1 600 

TAD I RWTAPE 

0008 

1236 

2200 

ISZ RWTAPE 

0009 

1207 

3327 

DCA RWBLK 

0010 

1210 

7040 

CMA 

001 1 

121 1 

1600 

TAD I RWTAPE 

00 1 2 

1212 

2200 

ISZ RWTAPE 

0013 

1213 

3326 

DC A RWCLOC 

0014 

1214 

1 323 

TAD EWM12 

001 5 

1215 

3322 

DCA RWTCNT 

00 1 6 

121 6 

1333 

RWBGN, TAD RWADBN 

0017 

121 7 

3731 

DCA I RWADCA 

00 18 

1 220 

1325 

TAD RWREGA 

00 1 9 

1221 

0341 

AND RW0077 

0020 

1222 

70 43 

CMA 

0321 

1 223 

3330 

DCA RWBCNT 

0022 

0023 

1224 

1 325 

TAD RWREGA 

0024 

1225 

0334 

AND RW7000 

0025 

1226 

1335 

TAD RWB610 

0026 

1227 

6766 

DTCA DTXA 

0027 

1230 

4330 

RWREV, JMS KUWAIT 

0028 

1231 

5240 

JMP RWFWD-4 

0029 

1 232 

1324 

TAD RWEsN 

0030 

1233 

70 40 

CMA 

0031 

1234 

1 344 

TAD Ml 

0032 

1235 

1327 

TAD RWBLK 

0033 

1236 

7710 

SPA CLA 

0034 

123 7 

5230 

JMP RWREV 

0035 

0036 

1240 

1325 

TAD RWREGA 

0037 

1241 

0334 

AND RW7000 

0038 

1 242 

1336 

TAD RW0219 

0339 

1243 

6766 

DTCA DTXA' 

0040 

1 2^4 

4300 

RWF WD, JMS RWWAIT 

00 41 

1 245 

5216 

JMP RWBGN 

00 42 

1246 

I 324 

TAD RW6N 

0643 

1247 

7041 

CIA 

3044 

12 53 

1 327 

TAD RWBLK 

0045 

1251 

7450 

SNA 

3346 

1252 

5256 

JMP ,+4 

004 7 

1253 

7710 

SPA CLA 

0048 

1254 

52 1 6 

JMP RWBGN 

■JO 49 

1255 

5244 

JMP RV.FWD 

00 53 

33 51 

1256 

1 325 

TAD RWREGA 

00 52 

1 257 

71 1 2 

CLL RTR 

30 53 

1260 

70 1 0 

RAR 


Standard RWTAPE routine per 
Jim Crapuchettes. Modified to 
give 2-block turn-around for the 
benefit of n^iddle-aged TU55s 
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AUTOCORRELATION ROUTINE MAY*1973 
RECORD 6 


03 54 

1261 

03.42 

00 55 

1262 

1 337 

0056 

1263 

6764 

03 57 

1 264 

1 326 

00 58 

1265 

3731 

0359 



00 60 

1266 

2330 

0061 

1 267 

5273 

00 62 

1270 

1340 

00 63 

1 271 

6764 

0064 

12 72 

5600 

0065 

1273 

1343 

0066 

1274 

3732 

0367 

1275 

4300 

0068 

1276 

5216 

0069' 

1277 

5266 

03 70 



03 71 

1300 

7402 

03 72 

1301 

6764 

30 73 

1302 

6771 

00 74 

1 303 

5302 

0075 

130 4 

6772 

0076 

130 5 

7700 

00 77 

1306 

5320 

00 78 

1307 

2322 

30 79 

1310 

5700 

0380 

1311 

6761 

0081 

1312 

0340 

0082 

1313 

7440 

0083 

1314 

6764 

0084 

1315 

67 72 

0385 

1316 

7432 

0086 

1317 

5316 

0087 

1323 

2300 

0088 

1321 

5730 

0089 



0093 



0091 

1322 

0003 

0092 

1323 

7766 

0093 

1324 

0030 

0094 

1325 

0000 

009 5 

1326 

3000 

0396 

1 327 

0030 

0097 

1330 

C0O0 

0398 

1331 

7755 

0099 

1332 

7754 

0130 

1333 

1 324 

0101 

1 334 

7000 

0132 

1 335 

0610 

0103 

1 336 

0210 

0104 

1337 

0110 

0105 

1 340 

02 £3 

0106 

1341 

00 7 7 

3107 

13 42 

00 70 


AND R W wO 70 
TAD RWO 113 
DTXA 

TAD RWCLOC 
DC A I R WA DC A 

RWAGN* ISZ RWBCNT 
JMP .+4 
TAD RW0200 
DTXA 

JMP I RWTAPE 
TAD RWM200 
DCA I RWADWC 
JMS RWWAIT 
JMP RWBGN 
JMP RWAGN 

RWWAIT* HLT 
DTXA 
DTSF 
JMP .-1 
DTRB 
SMA CLA 
JMP RWWOUT 
ISZ RWTCNT 
JMP I RWWAIT 
DTRA 

AND RW0200 

SZA 

DTXA 

DTRB 

HLT 

JMP .-1 

RWWOUT*ISZ RWWAIT 
JMP I RWWAIT 


RWTCNT*0000 
RWM12, -12 
RWBN,0000 
RWR EGA *00-30 
RWCLOC*0030 
RVJBLK* 0303 
RWBCNT*3330 
RWADCA*7755 
RWADWC*7754 
RWADBN*RWBN 
RW7000*7330 
R W 3 6 1 9 * 0 6 1 3 
1< WOP 1 0* 02 l 0 
RW0110*0113 
F.WO203* 3230 
R WO0 7 7 * 0 0 7 7 
RW0073 * 03 73 


T. GLATTKE 
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AUTOCORKELATIOi'i ROUTINE KAY,1973 T. GLATTKE 

RECORD 6 


0108 1343 7600 RWM200,-200 

0109 1344 7777 Ml, 7777 

01 10 
01 1 1 



f 
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A! ITOCORRELATT ON ROUT I NE 
RECORD 7 


MAY,1973 


T. GLATTKE 


000 1 
000!? 

1 400 

7000 

* 1 400 

XYPROD, NOP 


3003 

140 1 

7300 

CL A 

CLL 


0004 

1 402 

3130 

DCA 

XY1 


0005 

1 403 

3131 

DCA 

XY2 


0006 

1 404 

3132 

DCA 

XY 3 


0007 

1 435 

3133 

DCA 

XY 4 


0008 

1 406 

4031 

JMS 

RESET 


3009 

0010 

1 40 7 

1020 

DOXY,TAD POINT 

Setup delay 

001 1 

1 4 1 0 

3 030 

TAD 

DELAY 


0012 

141 1 

3326 

DCA 

DPOINT 


0013 

1412 

1926 

TAD 

DPOINT 


001 4 

1413 

10£7 

TAD 

CHECK 


0015 
001 6 

1414 

1415 

7510 

5222 

SPA 

JMP 

. +5 


031 7 

1416 

7390 

CLA 

CLL 


0018 

1 41 7 

1026 

TAD 

DPOINT 


0019 

1420 

1 023 

TAD 

RCOUNT 


0020 

1 421 

3026 

DCA 

DPOINT 


0021 

1 422 

7200 

CLA 

CLL 


0022 

1 423 

1 423 

TAD 

I POINT 

Put contents of location indicated 

0023 

1 424 

31 10 

DCA 

XI 

by "point” into multiplicand 

0024 

1425 

3111 

DCA 

X2 

3325 

1 426 

31 12 

DCA 

X3 


0026 

1 427 

31 1 3 

DCA 

X4 


0027 

1 430 

1 426 

TAD 

I DPOINT 


0028 

1 431 

3114 

DCA 

VI 

And the "delay” data into the 

3029 

1 432 

3115 

DCA 

Y2 

multiplier 

3033 

1 433 

31 1 6 

DCA 

Y 3 

0031 

1 434 

3117 

DCA 

Y 4 


0032 

0033 

l 435 

4536 

JMS 

I MUL 


03 3 A 

1 436 

7200 

CLA 

CLL 


0035 

1 437 

1 104 

TAD 

0R1 


003 6 

1 440 

1 1 30 

TAD 

XY 1 


0037 

1 441 

3130 

DCA 

XY 1 

And accumulate the XY products 

30 3 S 

1 442 

7004 

RAL 

. 

0339 

1 443 

1131 

TAD 

XY2 

in the XY registers 

0040 

1 444 

3131 

DCA 

XY2 


03 41 
0042 

1 445 

1 446 

7094 

1 1 32 

RAL 

TAD 

XY3 


00 43 

1 447 

3132 

DCA 

XY3 


0344 
004 5 

1 450 

1 451 

7904 

1 1 33 

RAL 

TAD 

XY 4 


30 46 

I 452 

3133 

DCA 

XY 4 


0047 

1 453 

7300 

CLA 

CLL 


00 48 

1 454 

1 1 0 5 

TAD 

OR 2 


0049 

1455 

1131 

TAD 

XY2 


30 53 

1 456 

3131 

DCA 

XY2 


3351 

3952 

1 457 

1 469 

79 9 4 

1 1 32 

RAL 

TAD 

XY3 


0353 

1 461 

3132 

DCA 

XY 3 
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A U T 0 C 0 h K £ L A T 1 0 N R 0 U TI \' E 
RECORD 7 


MAY t 19 73 


T. GLATTKE 


0054 

1 462 

709 4 

RAL 


0O55 

1 463 

1 1 33 

TAD 

XY4 

0056 

1 464 

31 33 

DCA 

XY 4 

0057 

1 465 

7399 

CLA 

CLL 

0958 

1 466 

1 106 

TAD 

0R3 

00 59 

1 467 

1 1 32 

TAD 

XY 3 

00 60 

1 479 

2132 

DCA 

XY3 

00 61 

1 471 

7094 

RAL 


00 62 

1 472 

1 1 33 

TAD 

XY 4 

09 63 

1 473 

3133 

DCA 

XY 4 

0364 

1 474 

7300 

CLA 

CLL 

3365 

1 475 

119 7 

TAD 

OR 4 

0066 

1 476 

1 1 33 

TAD 

XY 4 

00 67 

1 477 

3133 

DCA 

XY4 

00 68 

1 509 

23 23 

ISZ 

POINT 

03 69 

1 531 

2322 

ISZ 

COUNT 

03 70 

1 502 

S29 7 

JMP 

DOXY 

03 71 
0072 

1 503 

7330 

CLA 

CLL 

0073 

1 534 

1 1 33 

TAD 

XY 1 

00 7 4 

1 53 5 

3110 

DCA 

XI 

037 5 

1 506 

1131 

TAD 

XY2 

00 76 

1 50 7 

3111 

DCA 

X2 

0377 

1510 

1 1 32 

TAD 

XY3 

03 78 

151 1 

3112 

DCA 

X3 

0379 

1512 

1 1 33 

TAD 

XY4 

0380 

l 513 

3113 

DCA 

X4 

058 1 

1 51 4 

1023 

TAD 

RCOUNT 

0382 

1515 

70 41 

CIA 


3383 

1516 

3114 

DCA 

Y 1 

0084 

1 517 

3115 

DCA 

Y2 

008 5 

1 520 

3116 

DCA 

Y3 

038 6 

1 521 

3117 

DCA 

Y 4 

0337 

1 522 

4536 

JMS 

I MUL 

0338 

1 523 

1 104 

TAD 

OR 1 

0089 

1 524 

3133 

DCA 

X Y 1 

0093 

1 525 

119 5 

TAD 

0R2 

0091 

1 526 

3131 

DCA 

XY2 

0092 

1 527 

1106 

TAD 

0R3 

0093 

1 539 

3132 

DCA 

XY3 

039 4 

1 531 

110 7 

TAD 

OR 4 

539 5 

1 532 

3133 

DCA 

XY 4 

9396 

1 533 

5693 

JMP 

I XYPROD 

9997 

3098 

l 600 

7300 

* 1 690 

DIVIO, NOP 

0 g i) 

1 691 

7300 

CLA 

CLL 

9103 

1 692 

1 1 39 

TAD 

XY 1 

9101 

1 693 

31 10 

DCA 

XI 

3 132 

1 60 4 

1131 

TAD 

XY2 

3193 

1 60 5 

3111 

DCA 

X2 

0 1 0 4 

1 69 6 

1 1 32 

TAD 

XY3 

•0 105 

1 63 7 

3112 

DCA 

X 3 

'3 1 0 6 

1610 

1 1 33 

TAD 

XY 4 

•31 3 7 

1611 

3113 

DCA 

X 4 


Done? 


Then obtain by putting 

the IKV into multiplicand 
and N into multiplier 


The dividend is formed as in the 
case of the divisor, but the SIGN 
is noted. 
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AUTOCORRELATION ROUTINE MAY,1573 T. GLATTKE 

RECORD 7 


3108 

1612 

1 1 00 

TAD 

SUM 1 

0109 

1613 

3114 

DCA 

Y 1 

31 10 

1614 

1 131 

TAD 

SUM2 

01 1 1 

161 5 

31 1 5 

DCA 

Y2 

0112 

1 616 

1 102 

TAD 

SUM 3 

0113 

1 61 7 

3116 

DCA 

Y3 

0114 

1 620 

1 103 

TAD 

SUM 4 

01 1 5 

1 621 

3117 

DCA 

Y 4 

011 6 
01 1 7 

1 622 

4541 

JMS 

I SUBTR 

01 18 

1 623 

1014 

TAD 

SIGN 

01 19 

1 624 

7650 

SNA 

CLA 

0120 

1 625 

5246 

JMP 

UNNML 

0121 

1626 

7300 

CLA 

CLL 

0122 

1 627 

3315 

DCA 

NEGATE 

0123 

1 630 

1 104 

TAD 

OR 1 

0124 

1 631 

3110 

DCA 

XI 

0125 

1 632 

1105 

TAD 

0R2 

0126 

1 633 

3111 

DCA 

X2 

0127 

1 634 

1 106 

TAD 

0R3 

0128 

1635 

3112 

DCA 

X3 

0129 

1 636 

1 1 07 

TAD 

OR 4 

0130 

1 637 

31 1 3 

DCA 

X4 

0131 

1 640 

10 16 

TAD 

K 1 900 

9132 

1 641 

3114 

DCA 

Y 1 

0133 

1 642 

31 1 5 

DCA 

Y 2 

0134 

1 643 

3116 

DCA 

Y3 

01 35 
0136 

1644 

31 1 7 

DCA 

Y 4 

0137 

0138 

1 645 

5303 

JMP 

XDO 

0139 
01 40 

1 646 

1 647 

7300 

7001 

UNNML# CLA CLL 
I AC 

3141 

1 650 

39 15 

DCA 

NEGATE 

0142 

0143 

1 651 
1652 

1 1 04 

7041 

TAD 

CIA 

OR 1 

01 44 
01 45 

1 653 

1 654 

31 10 
7004 

DCA 

RAL 

X! 

01 46 

l 655 

330 5 

DCA 

LINC 

0147 
01 43 

1 656 
165 7 

1 105 

70 40 

TAD 

CMA 

0R2 

01 49 

1 669 

1395 

TAD 

LINC 

0153 

91 51 

1661 

1662 

3111 
730 4 

DCA 

RAL 

X2 

0152 

1 663 

339 5 

DCA 

LINC 

01 53 
0154 

1 664 

1 665 

1106 

7349 

TAD 

CMA 

0R3 

01 55 

1 6 6 6 

1305 

TAD 

LINC 

0156 

M1 5 7 

1 667 

1 6 70 

3112 

7904 

DCA 

RAL 

X3 

3 1 L i 

1671 

33-3 5 

DCA I 

L I NC 

01 59 
0160 

1 6 72 

1 6 73 

1107 

79 43 

TAD OR4 

CMA 

0 1 6 1 

1 6 74 

130 5 

TAD LINC 


If signal, the resulting r -will be 

positive (or 0). And the denominator 
is simply multiplied by lOOOg to 
provide for suitable display of the 
results. 


If the sign is 0 , the result 
'will be (™) and the correlation 
is obtained by multiplying the 
2s comp of the dividend by lOOOg 
...as above... and keeping track of 
the (-) indicator. 

r 
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AUTOC 0 RRELATION ROUTINE 
RECORD 7 


MAY » 19 73 


T. 6LATTKE 


9162 

1 675 

3113 

DC A 

X 4 

0163 

1676 

1016 

TAD 

K 1 000 

0164 

1 6 77 

31 1 4 

DCA 

Y1 

0165 

1 700 

3115 

DCA 

Y2 

0166 

1 701 

3116 

DCA 

Y3 

01 67 

1 702 

31 1 7 

DCA 

Y 4 

0168 

1 703 

4536 

XDQ* 

JMS I MUL 

01 69 

1 704 

5600 

JMP 

I DIVID 

01 70 

1 70 5 

0000 

LINC> 0000 









( 
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AUTOCQRRELATI ON ROUTINE 
RECORD 8 


MAY »19 73 


T. GLATTKE 


003 1 

BEGRLK 

0032 

BLKS 

0903 

CHECK 

3004 

CORCNT 

0035 

CORE 

0006 

CORK 

9007 

COUNT 

0908 

CROSPD 

0009 

DELAY 

0019 

DENOM 

00! 1 

DENOM1 

0012 

DEN0M2 

0013 

DENOM3 

00 3 4 

DENOM4 

001 5 

DISPLY 

0016 

DIVIO 

0017 

DIVIS 

0015 

DO 

0019 

DOCO 

0020 

DOXY 

0021 

DPOINT 

0022 

DXNT 

0023 

GET 

0024 

GETX 

0025 

GETX2 

0026 

HALF 

0027 

INS 

0028 

I NCR 

0029 

KEEP 

0030 

KPOINT 

0031 

K1 000 

0032 

K2 

0033 

K 4 

0034 

K 7 

003 5 

Kill 

0036 

LIMIT 

003? 

LING 

0038 

MINUS 

0039 

MPY 

0040 

MIJL 

0041 

Ml 

0042 

NEGATE 

09 43 

NEXT 

0044 

NUMER 

03 45 

OR 1 

0046 

0R2 

0047 

0R3 

0948 

OR 4 

00 4 9 

OUT 

00 59 

PN TR 

0051 

POINT 

39 52 

RCQFK 

03 53 

RC 0 UN I 


t 


09 1 0 

30! I 

0327 

0013 

00 46 

0324 

G022 

3263 

3033 

01 40 

0124 

0125 

0126 

0127 

0400 

1 600 

1 122 

0633 

0270 

1 497 

0026 

0 533 

04 72 

0 534 

1030 

00 79 

0477 

0146 

0152 

01 50 

0016 

0162 

0532 

0161 

001 7 

0061 

1 70 5 

0 7 1 6 

3 600 

3136 

1 344 

3315 

0351 

3143 

0134 

0 1 9 5 

0 1 0 6 

0137 

6454 

1016 

30 29 

0025 

0923 
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AUTOCORRELATION ROUTINE MAY*1973 
RECORD 8 


00 5 4 

READ 

0134 

0055 

RESET 

0031 

03 56 

RESETC 

0517 

0057 

RESETD 

9 53 7 

0053 

RPOINT 

0021 

60 59 

RW 

0147 

03 60 

R'WADBN 

1 333 

0061 

RWADCA 

1331 

00 62 

RWADWC 

1332 

0063 

RWAGN 

1266 

meA 

RWBCNT 

1330 

3065 

RWB6N 

121 6 

3366 

RWBLK 

1 327 

3067 

RWBN 

1324 

0063 

RWCLOC 

1 326 

0069 

RWF W D 

1244 

00 70 

RWM1 2 

1323 

00 71 

RWM200 

1343 

00 72 

RWREGA 

1325 

0073 

RWREV 

1230 

0074 

RWTAPE 

1200 

3075 

RWTCNT 

1322 

0076 

RWWAIT 

1 300 

0377 

RWWOUT 

1320 

mis 

RW00 70 

1342 

0079 

RW0077 

1341 

0080 

RWO110 

1337 

0681 

RW0200 

1340 

0382 

RWO2 3 0 

1336 

0683 

RWO610 

1335 

0084 

RW7006 

1 334 

0085 

SETCOR 

0041 

0086 

SIGN 

0014 

0087 

SMX21 

0120 

0088 

SMX22 

3121 

3089 

SMX23 

0122 

0390 

SMX24 

0123 

009 1 

START 

0200 

0092 

SUBTR 

0141 

0393 

SUMX 

0135 

0094 

SUMX2 

0137 

009 5 

SUM 1 

0100 

0396 

SUM2 

U 1 01 

3097 

SUM 3 

0102 

0098 

SUM 4 

0103 

3099 

TALLY 

0 1 44 

9190 

TFLAG 

01 45 

0101 

UNNML 

1 646 

3132 

WBLK 

00 1 2 

3133 

XAXIS 

0151 

0 1 3 4 

XDO 

1 733 

9 1 3 5 

XPKOD 

0 1 42 

3106 

XSET 

01 53 

0107 

XSTART 

0 4 70 


T. GLATTKE 
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AUTOCORRELATION' ROUTIN'? 
RECORD 8 


MAY,1973 


T. GLATTKE 


t 


XYPROD 

1 400 

XY1 

0130 

XY 2 

0131 

XY3 . 

0132 

XY4 

0133 

XI 

0110 

X2 

0111 

X3 

0112 

X 4 

0113 

Y1 

0114 

Y2 

0115 

Y3 

0116 

Y 4 

0117 


0108 

0109 

01 10 

01 1 1 

3112 

01 1 3 

01 1 A 

01 1 5 

01 1 6 

01 1 7 

01 18 

01 1 9 

0120 
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